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A. Problems 


Thermal data acquired aboard a U2 NASA aircraft over Walker Lane, 

Nevada were forwarded to us for examination. The data illustrated 
several noise problems; coherent system noise, scan line dropouts, 

"salt and pepper" noise, and gyro noise. At the HET meeting last 
January, Dr. Barnes indicated these problems were being addressed. 

Hopefully, these noise problems will be corrected before the first data 
flights in April over our test sites. 

One other problem was presented at the HET meeting in Phoenix. 

That problem was that ground calibration data from White Sands Missile 
Range might not be acquired for the HCMM Mission to calibrate the 
0.5 to 1.1 and 10.5 to 12.5 micrometer radiometer channels. We have 
addressed the impact on our experiment if no calibration data are 
acquired. (See attached letter.) We would like to reiterate that lack 
of calibration will affect our absolute thermal- inertia determinations 
and hence comparison between sites and over the seasons. 

B. Accomplishments 

The prelaunch investigation has focused on interpretation techniques 
development. We have worked on two of the principal tasks. 

The first of these is further evaluation of our thermal model. 

Last year we investigated the magnitude of the errors incurred in 
producing a thermal-inertia map and expressed these errors in terms of 
uncertainty in thermal inertia. These uncertainties were then placed in 
geologic context by noting the separation in thermal-inertia values for 
various geologic materials (see attached reprint). We have now developed 
a new method of calculating the Jaeger Laplace Transform solution for. 
surface temperature which is more efficient on the computer than the old 
method. The accuracy of this solution and of two other solutions - the 
linear Fourier series and the finite difference - were then recomputed. 

This time the analysis was performed using 48 time increments over a diurnal 
cycle instead of the 10 increments used previously. The results are that 
the Jaeger solution is essentially exact, and the relative accuracy of the 
finite difference algorithm is approximately 5% and the linear Fourier 
is approximately 5%. We are in the process of extending this analysis of 
the three algorithms to determining the accuracy of each in calculating 
geothermal heat flux. The benefits of this type of study are twofold: 
one is that the geologic application determines the choice of algorithms, 
and second is that the study is indicating different possible approaches 
to determine the geothermal heat flux. 

The next area of study was the development of a procedure for 
performing geometric registration. This procedure requires reference 
points to be identified concurrently on a control image and on a "distorted" 
image. These control points are used to construct a geometric correction 
grid. The correction vectors for the grid are determined by a weighting 
function of the four nearest control points using an inverse distance rule. 
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Once the grid is generated, the geometric correction is applied to the 
"distorted" image using a bilinear interpolation routine. The radiance 
values are assigned to the new grid using the nearest neighbor method. 
These routines have been programmed into our computer and are operational 

We now have certain data products available to proceed with some of 
the other phases of the investigation. We have solar and meteorological 
data acquired simultaneously at two separate sites. These data will be 
analyzed to determine the transient effects on a local versus a regional 
scale. We also have solar radiance and sky radiance data acquired on 
the same day at altitudes varying from sea level to 14,000 feet. These 
data will be analyzed to develop a model of the variation of solar 
(direct and diffuse) and sky radiance with altitude. Another product 
we have acquired is a map of the southern part of the Powder River 
Basin showing sand/shale ratios measured from hundreds of exploration 
drill-hole geophysical logs. The measurements lump the top 65 
meters of geologic section, but this much of the sequence is exposed 
across the area in the rolling hilly surface. The facies shown by 
sand/shale ratios in the Wasatch Formation match the previously 
determined patterns of vegetation cover displayed in specially 
processed spectral-ratio images from Landsat. The existence of 
facies which had not been mappable in normal geologic field work is 
important information which will bear heavily on the interpretation 
effort we proposed for HCMM data. The lateral change across the 
study area, from dominantly sandstone to dominantly shale, should be 
visible in thermal data at 500 m resolution. The measured sand/shale 
ratios will be used as control data for calculating such ratios from 
gamma-ray survey data ^available at course flight-line spacing for 
the entire Basin. These ratios will be used in a later study to 
correlate with the thermal-inertia images. 

C. Significant Results 

The significant results from this reporting period are: 

Development of a new more efficient algorithm for calculating 
surface temperature. 

Determined that this algorithm is essentially exact and that 
the relative accuracies in determining thermal inertia of the 
finite difference and the linear Fourier series algorithms are 
approximately 5% for both. 

Developed the procedure for performing geometric registration. 
For details of these results, see the accomplishments section. 
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D. Publications 


Presentations were made during this reporting period by four members 
of our HCMM investigation team. All of the lectures were presented at 
the SEG remote sensing school for exploration geophysicists held 
January 23-27, 1978, in Denver, Colorado. The titles of the lectures were: 

"Remote Sensing Overview," T. W. Offield 

"Digital Enhancement of Geologic Information in Multispectral 
Images," D. L. Sawatzky. 

"Case History - Thermal IR," K. Watson 
"Detailed Techniques - Thermal IR," S. H. Miller. 

E. Recommendations 


We recommend that NASA acquire the radiometric calibration for the 
HCMM investigators. See the problems section. 

F- intentionally 

G. Data Utility . - 

No U2 or satellite data have been available to us. 
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KVAUIATIlV: C'K Al.HiMiiTMM.' KtMi r.XOl.Oi; U~AI. 7l!K!iM,M.- ItlKKTlA MAI'l'iriC 
. M. Hllb l- :iiiJ Ki'iiru-' li ;'ii 

II.;!. li.-.i U’li ! <••1 ! ;!ut‘v<'.v, I'l'iiV’-'f, i’i> 1 oi'rtilo 
A!!;!TRACT 

'■'he .'I'l’or:-. 1 1 ti,: ;i I - 1 1 ;i m.ip .ti-o oI’ Uii'.'l’ r,o:i'?f;il 

tvip':!; n'.ar.iii’inH'M! , :i:ia 1 y:; I , and nu'di-1 a I r!|' I I I'l mil 1 oi! . Ta nni-ii;!:'. 1 CPo- 

!.;a I :a- 1 nvaiK.v o:' t'rri>rr. , w.' .•xnia';:;: Uii'.’n In tt-riaa n|' ii;ici-rl;a I In 

l.lan’ina 1 Innrtla and mirTn-ai'*.* t ia.nto unnfrla I nl 1 *■;*. with thn 1 - I noft l:i valima 

ai' ,-t-aK. rU; raa tan- 1 a ! a . ".'hua Un.' ai'p.) Icat.lona and jn-aclleal 1 I ml t at, 1 ona of V.liO 
^.^■chn 1 1 1 i;a t. ra t.nd . 

/'I! .nn-oVa ai-n o;i 1 nn 1 at. nd ii.sliw: 1 hn panana.-t.-i- valiins api'pnpn ] a to to .a alte 
at tin; Hal't Rlvon, Id. AltlKMi,-,!i tli.-an in-nof vaUioa aoi-vo to llluatrati.' tho inad- 
nltiah.-a l.liat o.an bo o.xpoct od I'ror t tio t,ht-..'o in-noi-al typoa of oni-opa , t^xt!■npo^a- 
tlon to othon alt»,‘.a ahould bo dono iiaiint par.irnti'P v.ajnos pant 1 cii I a n to tho 
•a rea . 

ooaaiiporai.nit opi-ona Inti-odm-od by nrn 1 1 1 apoo t I’a I acaimlnp: ayatoma con-.monly 
p.inro I'i'oia a :io 1 ao-mpi 1 v.a 1 mit-t omporal ufo d 1 iToponee (MKAT) of O.IK foi' a I r- 
oraft. ayatoma t.o IK foi' a.at.elllt,.- ayatiaaa. I'lio poaii U 1 nir, nneorta 1 nt loa In 
tbepri.i! ! .nonl la I'an^ro rr.iin 15 "'ID ( tbori::a 1 -Inonl la unit)* to 15d Till. 

n'hroo auifayy lenipor.atLiia.- alporlthma wopo o va 1 nat.-.l : llnoai- Koiiidop 

aoploa, flnltt; .1 ! |■^•.,■p-■noo , and I..i|'laoo tp.inafopn. In teema ol' poaii 1 1 hit; oppopc 
I n tbopmal lno.”tI;i, the I.apl.aoi.'-t.panafopn n'.othod la the moat atnnipat? (I'bO TIU), 
tlio fopwapd f In 1 t'.'~d 1 f foponoo ir.otlioil la I ntopraod l.ate ( 300 TIU), .and the linear 
Foiinlop aopIoa mc-ttiod t ho loaat. .aocupate (h60 TIU). lIpKcyop, the two more 
(.■xaet liu.’ttioda poiiulpo mope ci.smpui.op tlm>.', and both lack the ability of the 
Koup.K'P-.si’P i oa a 1 iiop 1 t.hm tis lllur.tp.ate the plsyolca) :i 1 >'n 1 f Icanoe of Individual 
topma. by oompaplnt’; t lu' I'ppopa with the panno' of thopiaal Inert la« of poolo.flc 
r.atepl.ala, It la poaalblo to aeloct tlio no.st oo.-.t/erfoot 1 vo alc.oplthra fop a 
pai’tloulap application. 

Model ;• I mp 1 1 f 1 fi.at 1 on epnoi-a penult from three i'.oupce;'. ; tp.ans'.lent effeet.s, 
topoitP.aphy , and aunface eeatlnp: effeota. For example, flux oi' 55 w,aLt.a,''rn' 

(ei]u 1 va lent t.(i a water ev.apop.ai I on rate of l.b mra/ilay) woul<t ppis..iui;e a thermal- 
Inertl.a oppijp of .about fOO TIU. If no toporraphlc eonrect Iona aia' made, a lO'’ 
.aout.hwestopn alope c;,au:jos an i;ppop In therm.al Inentl;i of about 350 TIU. A 
homatllio aurfaeo co.atlnp one centlnis.'ter tlilck on a rock will ppisduee an oppop 
of ap'ppox I na tel y TOO TIU, whereas .a one-iiil 1 1 Imeter layer thleknoaa will have a 
neyllftible effect (.35 TIU). 



The total oy.'.t.tm eppona In thermal Inertia are placed In fjeolofjlc context 
by notlnc the neparatlon In tlu.'pni.al -I nertla valuen fop varlcua feol oplo mate- 
rials. For example, the therni.i 1 -Inert la oijparat loti between limestones and 
dolomites Is typically 1300 TTU. The error an.ilysis technique Indleater. that 
therma 1 - 1 nert la d I scr !ml n.at. Ion .should be possible bvttwoen black shall’ and (tabbi'O 


(separation of ap'pi-ox Im.iLo ly 'I50 TTU); discrimination between these two mate- 
rials war. not possible iislnp bands.at spoctral-i’eriectance data. 

■ F.rror.s Ir; therm.al Inertia can be translated into crn'oi'S In bulk der.sltj’ 
and moisture content. Therma 1- 1 nei-t la mapplnip from ali’ci’aft (NKAT = O.IK) has 
.'in accui'acy of about 1' In bulk density or oo’.il v.alent l.v , a sensitivity csf 0.5J 
In water eontent . From satellite (NKAT = IK) these accuracies ;iro 9' and 35 
respeeUlyoly . 

A pr'aetical evaluation of the error an.'ily3l.s l.s demon.slrated foi’ .alrer.aft 
d;it a acqu I red at R:ift Ulver, Td . We c.annot discriminate tl.e tuff from the al- 
luvium, h.ivlnp a lower thermal Inei'tla, ap.d the iava, haylnp a hif.her tliermal 

• TIU = 1 watt • .see'v'm"/K 
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lijci'l. 1.1 , without .1 t opof'.i’.iphlo uofivict loii, but wo c;in dl :ii‘ i' Itnl iiu t e the lav.'i 
riowa rrum t lie ulluvliira. 

1 . INTKOlUIOrrON 

Therma ! - 1 iieet I a mapi'Ii'lt la .a reoeiit ly developed t<’elinliiue th.it ;i 1 lew:; dl:;- 
ee t ir. I n;U. I ell ainetu: peeleple material:; 111, |2|, |31. I*). I.'>1. I h I . 'I'he piirpe:;e 
of t.lil:; paper 1:; to dl:;cu:;:i t he types and tiiap.n 1 tiules of err.'r:; eiieeiuiti'i'i'd In 
applying, thi:; terhnl.nie an.l to expre:;:; these error:; In a ite'olop. lea 1 ly meanlnp.- 
I'li 1 way. The error:; are of three major type:;: mear-uremi’iit error;;, .aii;. 1 y :; I a 

error:;, and model :; I mp 1 1 f I rat I on and a;;:;iimpt K'li errors. The m.ajor foriis of 
this t'.aper Is to develop a pror.edure to seleet. an alt'orltimi which Is approprl- 
.ate to the peolople app 1 1 e.a t I on and which minimises eompiiter time. 

In .addltli'ii tvi present 1 lip. the mat hemat 1 ea 1 basis for the three a 1 por I thm:; , 
wo provide a luimi'rleal evaluation with respect, to an exact solution. Kst Iroat.es 
of the errors due to measurement and model :; 1 mpl I f lc;i t Ion .are comhineil with the 
.an.a lys ! ;; errors, .ami .all errors are expressed .as uneerla 1 nt 1 es In t lierm.i 1 
Inertia. These oomhlned errors are comp.ared to the t hernia 1 - 1 no rC la v.alues 
determlnovi for some selected t',ool<’|.'.le sltu.atlons. As .an 1 1 lust r.a t Kin of the 
tin.alysls teehnl.iiie, the l-'our 1 ei— ser I es ;ilporlthm Is applied to t.he construct Ion 
of ;i t. lierma 1 - 1 iiert 1 a m.ap of an .area In the Haft IMver, Id. 

.?. MATHKMATICAI. ltA;:i."t 

Thermtil Inertia Is .a physle.al property of peolople materials, and Its value 
can htt di't. t'rml tii.vl l.'V meamirements of t.he tlliirnal varl.atlon of* surface ti.'inper.a- 
ture and tin; albedo of tin- surface. A physle.al model for the siirf.ace teniper.a- 
ture variation Is therefore rei|iilrcd. 

We assume one-.I I mens I on:i 1 , periodic he.atlnp of. a uniform half space of 
const. :int thermal properties. The (rround temper.it uro obeys the diffusion 
equ.atloti:. 



.a." V ( X , 1. 1 , a V ( X , t ) 

“■ 3X- .at 

where v(x,t) Is the proiind temper.at.ure at a distance x below tin' surface .and 
a time t, .and e Is the t.lierm.al d 1 f I'lis 1 vl ty . 

To solve this edii.tllon for the surf.aoo trmporaturo , we require :i form of 
t.ho solution appropr I .a ti' to the diffusion oquallon and the houiid.ary oondll;K'ii, 
which assumes cimserv.at Ion of etierpy ovor the diurnal cycle. l<«;fore loeklnp at 
possible forms of tin: solution, we in;ed to luvestlpato the boundary condition. 

If wr asstinio a boat -ba ’..anco cotulltli'n exists at the Siirf.ace, then the 
Incomlni:; heat fliixes--the Incident solar r.adlatlon (1), the sky radiation (.S), 
and the poothermal 1 ;;.at flux (Q) — must balance the outpo I up fluxes — the con- 
vective flux (0), the ovapor;itlve flux (B). and the pround reradlatlon flii.x (fi). 
The resulting expression of the heat flux «t the surface Is: 


■'I + S + Q-G-K-R {?.) 

x*0 

where K Is the conduct I vl ty ol' the m.ati'rl.il . Tho eonveotivo ;iiid ev.aporatlve 
flux;’S are Included In e.piat.len J for co;np!eteness'; howovov those terms will he 
handled only qualitatively In. this p.iper. 

The p.ai'anieter of Interest, the siirf.ace toniper.atiire , Is present In equation 
.a In the In’.at-coiuluct Ion term (the left hand side of equation P) and In the 
groiind-rer.ad I a t len ti.'rni (K). Kqu.at Ion P Is rewritten to dlspl.ay the surface 
temperatui'c : 



( 1 ) 
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3 V ( X , r. ) 


“ -cd v(0,l)‘* + a(t) ♦ Q 


w!»-ri- c I:-, tin- mt'aii om ! r.:-. 1 v 1 t.y of t ho j^foiiiul, o 1:-, Uio Stor.-iii- 
::tant, a ( t: ) I:; a t.t'rm comb I ri 1 tij': Ku' I. Imo-ilopotiilonl. aoiai’ CIiix , 
riu.\ and t.luia Inoluilo;' t li.' aurt'aco aU>odo and om^ sr. \ v 1 ty , and 
l.hormal boat I'l nx. Thia i-xpi’oaa Ion miiat. bo aat. lai'lod Vov any 
of v(x,t.) clioaon foi- t.ho aolullon l.o Mio dll'fiialon o.nialli>a. 
form of thla boundany condition provonts a illooct analyt ical a 
dtffuaton od<nit.ton. An exact immootcal aolutton can Ih‘ calcul 
the oomploxl ty , the dlfflcailty of 1 1 luat oat Init the phyatiial at 
tooma, and tlio compiitco coata all poovtdo Incontlvc for Invoat 
exact alpoolthma. 


Holtamann con- 
and aky radiation 
tj in the ff.eo- 
i.'ipiat Iona 1 form 
'I'hr nonl 1 near 
oliit. liMi of the 
atrd (71; however, 
p.ntftcanco of 
I gating loaa 


One method of handl<ng the nonlinear boundary condition Is to Ilnearlco 
It and to express the ac.utton to this condition In an exact form as an Infi- 
nite Fourier aeries. The 1 Inear 1 ;-..at Ion la .achieved by performing a Tay lor-ser I es 
expansion of the surface terap<M‘ature around the sky temperature and discarding 
quadratic and higher order terns. This solution has been discussed In detail 
by Watson (I|,,nnd the form of Mils solution Is: 


' '^sky * 0-1^) V 


t'os(not - 

n‘^0 ^(s+r^)? 4 (r^) 


where 3 

" = “-'^sky 


*n “ 


\ s + r vTT / 


and Is the sky temperature,. A Js the grovitid albedo, C is an effective 

atmosphei'l c-transmt ss Ion factor, Is the solai* constant , Q Is the gi'othcrm.al 

heat flux, and are the anplluudc .and phase of the nt_h harmonic of the 

Inso-latlon, and I’ Is the thermal Inertia. 

Another method of solving the diffusion equation Is to use the exact, non- 
linear boundary condition and to apply a flnlte-dl fference Iterative technique 
(91 . This solutloti assumes the form: • 

V - 2*V ♦ V = (V - V 1 (Si 

mtl,n m,n m-1 ,n At- k ' m,n+l ni,n' 

where V Is the ground temperature at a time, nAt, and distance from the sur- 
m,n ' ' 

face, mAx; Ax and At are small Increments In distance and time; and k Is the 
dlffuslvlty. 

Theforward f Inltc-d I rfei-enoe method requires an Initial solution and pro- 
pagates thl.s solution foi-wai-d In time. To Insure eontlnulty at the boundary, 
the heat-balance boundary condition must be satisfied. The boundary condition 
for the forward finite-difference solution Is expressed as: 


^.n.l " ^ '■{''0,n+l>(-'*l’ 


where K Is the conductivity and f Is the flux (cq. 2). 
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ii" 1 . . . . , 


' ‘n '■■ ' ilii.c !:!(,• tin- fii-oiii’a I :i i li,- nth t I i:!i- 1 tit i-fta 1 ; 1' 1 a tlio 

. la ; h- i-t-l.M .,r M,.- hfatli,,: rh:x; v. la . h,. a...T;..-o 

iMT'ni'V”;''. ’’’ V"’.”- <1^- t.:-,. of ,nt„„.flo:tl ooof- 

h-^- ,/l r".i'V *'•' ''!• Itu.-rvala In A ,1,-! a J 1 o.l' 

o'-"', lv,f'-','^' •'•■•■t-o.- |71 for ll,.. 0 ,-,-UvlIo honlltn: 
layi^- otioi/'a half al,!-;. ' "■'■to.i. by K:u a.-tt (SI to fho hoaUtht of a’ 

A. AMAl.Yf!;: Kh'KOH;: 

I'foaonl a Uio iita Uiomat I o.a 1 ha::!.:; fo,- t.ho t hfoo al.-o- 

''•b'j'o iiao.l to conirnto tho.dln.-nal furfaoo I .auporat uro.-. jt ) l„ thin 

-;o. .ton, Ko nuiiiofloa! 1,Y ootr.t-afo i ho orfora hy oxi'i-oa.-.ih,.: : ho t ."iirofauaa- 1 f- 
I, i.’tiooii a.-, oiiulyalont ,il ffot-o, S„ t horir.a 1- Itioihla. fho .1 1 f f 1 ou 1 1 v In r,'f- 

00. \ tift 1,0 phyaloal ,-• ! ,:;i I f 1 caa.to of 1 ofna ami t ho amotmt. of o.'iiiput i lir.’ 

fO.Hitfod for oaol, alitoflt.hr,, at-o alao .llaouaaoA. ' ui . . 

Tho ovalaat Ion of u.o a I, top I Uinta In , ot-nia of nn.-opt.a 1 „t 1 oa In tho-nnj 
nopt. a a poU-oPmo.! In t wo afopa. Tho final la to' non, nano fh;. fhnoo a 'onllhtna 
wtth oa. h Ot.hon by appnoxlntaUf.f. t.ho luoMonl aola,. aa a atmiaot.tal haU' 

• wavo. ,Tho a.vott,! atop la t;n .-ompap.' t.ho moat, oxaol a 1 ,n-,M t hlLlfho 1 ■ 1 1 . no 

1. pana.t.' Pin— wit h an oxaot t h.-op.-l 1 pa 1 aoVni.ton— a pnno alnna,.Ul flnv. ri,l.; . 
oatabt ah.ta t-lu> PolaUonahlpa of thc llu-.'.- tnolhoja t.o U,.- oxant aolnt lon.' The 


Laplant'-I. pana form A I for I l;lim. 
V! n 1 t.;-il I f fop..,n;,t A lftoivlviira 
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Till 
31H7 Till 
'ifiO Tin 


Table 1 plvea Vh.t Hat of pap:, note, -valuea naod to mat.; tlw above .tonipanl aon. 

Tlut polatlv.’ .'p.lopIn,< .of Ui.f a |, top ! t. lima wlti, p.?a|-ct to anotiPinv la i- 

alran-iu.m, whleh la an ImpHnlt tv lat I onat, I p 
hetw.'on M.pla.to t .-mpopa t. npo ami aupfao..- fl„x, la the moat annupato- noxt tho 

p.-al'M-" ^'p^ ■""*• '•'"■"•J'-'-^-->-'-lea a 1 .to,' 1 1 Im, , wlVloli 

111.. .11 1... .. u,.y I .I.il.,tl.tri-I Itmai.-p t.-i-ma. The n,,a,;n 1 V lule of M,.' .-ppopa nan ho 
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An Important nona lOora 1 1 on In tho eoat./efr.-ctlvoneaa of t hepm, 1 -1 nopt 1 a 
mapping la , h-t .t.vnpnto,. time P...,utp.M te aolve the vaelona :'lg.l-irhma; ll U,e 
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:ir..» i nr ro\nt K't'n-'oKr:* U'.t to 


Wat on (I) .iovolortnp mot ItOvir t'o-r Iv* h t !n'r:r.:i ! ! r.oin I .a-::inpp I n** ;i::.l 't hi'rm.a I 

h*vit-tMiJ\ mnrr'no. 


KaMo :o;mm.ar K.-.o;- t h.«* .-I'niltj; ot* thin roottv'n. 

Tv'-’VAK .''VtVVKN KHtU'tiit 

y.i'.arnt'omor.t ■ error;', t nt ro.tuot'.! t'y •!*.<' t:u: : t ! :*-root r.al :*o;ir,n!r»: r-yo-t**::; l::-ro;'o 

.a t*{'.r..Kv.*.'.or;t n 1 Ir.:-*, rnr.t'r,*. n 1 I K;;'. t at *. on ot; t !:or::;n t - nt-i't Ki r.npp K*.: . . '■■lio;-.o orr.'r:'* 

.arr t h.i* roaoit o>'::;h tr-.it .'f :'.K*:or;*, t r..* t t h.o lU'Ku' l!i!5tt of t ho 

.iotootor, .'i;‘;r : ! o!*r , t.aro rooor.'or, n r-nt o ! t ! : *' ;-.v:'.t t'r , ii*’wn- i ! ::k 

t o I ortot 'oy . ‘.’tn'j'o . <'rror;'« :ir*' .'.'rr'^oo t y «'\pro;von o,;. r.o t - m * v:> *; or.t - 1 «'mporat tiro 

si ! t't'oror.oo;’. K.'r ;:tr,--M:‘* :*-y:'.t *•::*, a it valtn* t*s'r n t a t o -.‘ 1 ho 

.at't ayr'i’ri ! , 1 K . A :K NKaT ! ;• vor*'o::tnit a t t *. • -'t' r.nt'olilto ;*i:oh mj** 

NKiihna i I K an.i ;V; the IK'MM oHoat .‘ara.-tty y.iroKu- :'at'o:':'to ryatom 

may ’‘..-wo a :'• ! ! p.r.t ] y Knvi'r '.ii-.U' v.i h:or. .'t* ;,•• 'ort a i rt y ' r. t oi-.rorat ::r«‘ 

s's'nvt'rt t *» .arpt',' X ' ;::a t o 1 y 5'* Tit' tk'f t ho I'.oorat**. .*oi.i 'VH't'or a iKitoIt't**. 


".‘hiO; 

»' 

r.'io.. 


a: 

o 

In .1 


1 . 

r. .1 ,• 

\ ho 

•o'-.T 

a thtr 

s'sh 

s'O.' 

I'.v 1 t!i- m 

irp I np 

i 

r. ' 

I'.n' 

;i;U r>'rr' 

*ao5il a 

t* 

un.l.'ini'tu 

•1 ! 

t i:i 

1 ! at t 

' n t : 

heron* 

* O 

! 

t ‘ t Os' 

I'.!!!.; 

lO . 


A 

n oi 

i' * 

'h t o r 

np 

Pli!,- 

:lat 

:\ 

with 

tn't h 

.1 t tv 

V.’lH M 

nil 

a.M *s 

olUt 

I- .•.!■. 

t ttat 

rr.'.i 

i:.‘ 

Oil ! 


*0, V* 

:* 

f'T.I 

.••t . r 

n 

O 

ii's';'. 

ts'.p.:'. 

!*o.l at 

*il 

lY.i 

ror.li 

! i5f 

*? . 

Ih'i'.atij’- 

o or 


In t-iio 

V.'lll.' 

to p 

at 

!i an 

1 .aoannor 

.'!• \ O'l* 

■i! I 

s';i 

(pi t r 

ii. y 

iw, nnit 

j*.'i 1 


1 

.'in! 

ho pr*' 

un 

.1 wl 

t \ n 


h..» % 

1 >■,.■ 

asimo 

I .ma pa' 

p ‘ 

ora 

‘ nat {. 

;i M ! 

tW.' 



1. •. \ 

ii llr.ii';!. 

I'ayt I 

mo 


ar 

nu\r 

.-4 a; 

1 }sj. 

aam<* i 

Imo 

at 

o nut. 


K'al ly ' 


:!! 

or 0*1 

■ h Uli i-;i. 


r 

ar.st 

a ro 

tv*’ 


i-l.-.i! 

ly Ki 

otn ls*a 

1 

"'!u 

ao "it 

;l!. i'l'.i'it" 


1 .* ;? 

.a Pi’ 

.•nlj.-.i 


at 

ra. 


• r 

-rii> 

! I'.'lt 1 

•n ia 

ro: 

pm* 

.1 I 

V :n‘i 

*•*'*. 1 

■r. 

r«T t*l.’ 

ro 

iU* o 


ot\ a t: 

:u« 

t or 

tma,: 

o» 

vnll- 

.•uul 

OsUao 

psontl i 

y .1 

s.ioi 

■sKryl 

np t 

is'ao 

!\‘.a t. 

ur 

or, on 1 hat ror 

\ n I no 

** 

i t :•*. 

s'l't o 

.1 ** 

1 ni.ii'.i'-t’.i 

Irr 


; h.« 

Ml 

pr.t 1 1! 

»' I h. 

*rm I 

Imap 

or 

. H 


I'oat u! 


tvio 

1 r 

! <• 

sf s'n 

hot h 

\ !u' 

rsi.ar', ot* 

rc 

V \ I 

s't ano 

*• Im 

ip.o an*l 

t ‘u- 

m:\ 

;-t I'r 

t-horma ’ 

Im.air.o 

V.Vsi 

at. h 

o K1 

or 

tin 

s-.t s>:i 

1 

’MUUs- 

r s 

sK' 

nl,-iu 



t h.oj 

ma 

i in 

a^-n'?*. init 

in.ay V** 

I 

osn-i! 

Osl s' 

:i 

o I 1 h 

sT t h 

• roi* 

lo.-t ar 

Or' 

. -r 

! ! : o t 

!ior:r. 

it Irnafn* 

» > :* •' 

*1! 

! s' : 

to.l” slay 

! mo 1 1 : 

•Kt 

o-pa 

tr;;. 


T!:o 

oj'.s'to 

■ sK p 



wl 

o t 1 i « ' 

s h. o 

appoar- 

.atu'o 

s' 

r t lu' :*oaia;rt 

■ ia ,-s- 

Vo 

riu'.l 


t s' 

V'op.r 

arh.y , 

t tu*r 

sual ’.r 

*'rt 

Kv, 

s'r ; 

of ; s' 

s' l ;ir.s's* s 

The;: 

*• 

r«* t*r 

rons’o :*o; 

n 

ar 

i' 1 h 

on u 

a o 

.! a a 

: ::pn 

* 1 s' 

.a trl.i 

:*p. 

: a : 

i o 

1 

n tor 

M 1 >:o 

iv! 

! lim . 

















^nr.* o 

rr.'j’a ,i;!i 

’ t . ' !:; 1 


s'pta 

1 rat 

I,' 

n ar 

*■ s'T 

W*' 's 

ypo;-.. 

■V! 

s* t 

' ! ■ ;*. t 

* ;5 a 

ml 

rorM 

a t, 

!-;li ! 

on ot‘ a:i 

Kna r.o- 

ra 

tr t 

s' I ho 

m.ar*t 

s'p tm 

poo-p 

:i Ir, 

Thl 

a i: 

I ;; rop 

! t 

(t, 

I'l-l'-ul ( 

I'.i 

an , 

l;u l.'ii 


vs'rr 

t's' 

all't'sK'a 

.tr.si' 

s'l- Ir.s' 

s'ri 

s's' t 

t orr. 

*' pa i 

.r.*t'~ 

.il IT 

t' r 

•ni'i’ 

s's'ir.put :t t 

K'ii--!- 

v't 

!i <' t* 

will 

-.Is 

I's*;; 

■ si! t 

: .an 


t 

O s 

a 1 .'li : 

a! K- 

l s't* 


t h.orsr.il tr.ori Ka, 





The second type of error 
Inei-tla to a topof’raphlc nap 
a UuM-mal inertia to a par 
ttie mapnltude of these error’s; Itowcver, 


is a tnlsref;l strat Ion of a calculated thermal 
...... This error Is an Inaccuracy In the asslp.nnient of 

tlcular p.iropi’aphlcal point. We liavc not estiiivated 
we have i]ua.l 1 tat 1 vel y cxaiiilned them. 


Mloi’cp.l stratlon of either type results fi’orn one of two conditions; the 
spatial frequency of the nilsregl st I'a t Ion Is higher than the spacing of the I’e- 
fei’ence features Or- tlie itaage features are not discernible enough to select 
suitable i-efei-encc points , The first of these condltlons--the rapidly varying 
scene — Is cori’ectable by choosing tnore frequent refoi’cnce points, and the re- 
sulting ei-rors In thermal ' Inert la can be minimised. Tito • add 1 1 lon.r 1 points do 
impose a significantly higher computer time. The second condition — low scene 
contrast — Is common especially In nighttime Imagery. In such situations the 
best estimate of the location of a reference feature oit an Image-pair is made 
by scaling tbe distance from other’ Identifiable features. In this case, the 
resulting roglstrat Ion may be geomet”lcally Inaccirrate but numerically cqulva- 
lerrt duo to the resulting low contr’ati: image, which will not c.xhlblt much vari- 
ation in temperature and/or r-eflectaric’.' valiics. Thus, these errors pr’oducert 
in computing or assigning thermal Irrert ar’c small. The edge effects that may 

be produced by identifying fcatur’es b.’ . Ing can be reduced by spatial fil- 

tering. 


Tn addition to the measuremerrt er’r’or’s, arrother’ Important source of error 
in thermal -inertia mapping results from simplifications Introduced by a par- 
ticular model. Those s 1 mpl 1 f lea t Ions Include the exclusion of transient fluxes, 
topography, and/or surface-coating effects. 


Transient effects are those resulting from such conditions as cloudiness, 
windiness, or ground moisture. Effects of this type can bo viewed as a flux 
pciturbatlon on the pei’lodlc solution 1111. To Illustrate this situation, we 
con.slder a heat-flux perturbation of 35 watts/m^. The magnitude of this per- 
turbation Is equal to 3.5? of the solar constant; or an equivalent change_ln 
sky temperature of lO'^K, or the effect of an evaporation raVo' of 1 . •! x 10~® 
m/sec (1.3 mm/day). The translciit-ef feet calculation was made by varying the 
duration of the perturbation and noting the resulting change in suPfacc tem- 
perature at a later time’. A graph expressing the resulting errors in thermal 
Inertia Is presouted in figure 1. The maxim' m error introduced by this per- 
turbing flux is approximately 300 TIU. 


Correctable topographic effects require an assumption as to the scattorlug 
law Tar the sloping sui’facc. and the need to neglect reradlatlon botwet i ad.Ia- 
■ ”>t .•''orfacos . We shall assume that the surfaces are Lambert reflectors, and 
tiras the local solar flux is proportional to the cosine of the angle between 
the surface normal and the solar radiation. In that case, an oa-st or west 
slope to the surf.icc causes a phase shift in th diurnal temperature curves, 

.and a north or south slope to the surf.ace pi’oduces a change In amplitude of the 
dlui’iial curves. A surface oriented in any oth’ar direction produces a combina- 
tion of phase shift and amplitude changes. Both changes result in erroneous 
day/nlght temperature differences. Similarly, slope and surface orientation 
affect the reflectivity and can cause erroneous determl nat long of albedo. Both 
temperature and albedo errors contribute to an Incorrect determination of 
thermal Inertia. The magnitude of this error (using the Raft R'ver parameter 
v.alues) for a lO** slope facing southwest Is approximately 350 TIU with the max- 
imum error in thermal inertia occurring for. a southwestern facing surface (fig- 
ure 2 ) . 

Surface coatings change both the. surface reflectivity and the diurnal tem- 
perature; hence they affect the apparent thermal Inertia of geologic materials. 
We examined the thei’mal effect of a hematite coating over a half space, having 
the average thermal inertia of Igneous rocks . Elgur’c 3 shows the implied 
change in equivalent thermal inertia for various thickne’sses of this layer [81, 
For example a hematite layer of 10 millimeters will produce a change in. thermal 
inertia of approximately 300. TIU. If the layer thickness is less than 1 milli- 
meter, the coating is transparent in terms of thermal-inertia mapping but will 
be observable as an albedo difference. 
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or Uu' U.f.-,- r.lnpl !n.-:n Ion oriVot.;-. a.n-.r :• IK-d . oiilv t !u- t .vo.-npt'v >-o>- 

C ion in ;,t ;i nt.i,-:.' ot vlrvolopnoiu. to In- inoorp‘'r:iv<\t Into t lio ihofWri - 
Inoitln ninpp 1 Ilf, . Tin' loo.i! v;i ;• 1 ;iM onn In t !io niii-rnoo nt ;i ! nJ n>: , vofOt at Iv'n 
oovon, anO initfaoo iMitflino:'.:-. ootnir at. a\:u-alo whK-h pi<v 1 lalt :; , total lo.i .-loJ.'Miif 
booain-.o thtit, KotiM i-o.ni!i-o KiioVum.-.o aP t hono ol'Iivt:; at ovoiv noint ' 

• nt [to tout nlto. V\o oonnKtot- tl-.ono t..v} n of otToots oooiivf I tij', at a looal 
.■>a!o l ffoa - 10 U' t o tof praot.Ioal tino of ronot o-i-oin- ! n>" toohnUpion. o\wl lao- 

. t opof.i-aphio oitoo’ ,; j,, [ .fa 1 I ;:o,i at- .•an itiio ;.o ilopartiiro i-. om l.atitl'i'ft: oiiitn- 
nlon, nof.i.i lat. Ion pot woon atl.Iaoont. niirfaoo.-; , or nlruirw I no, offool.n aro a 1 no 
pi'oPaPl.v In t.lio naino oati-f.or.v. 

P. i'lPt'KlMIN'.ATl.ON Of-’ llK.’LOtllC KATKK I Al.t: 

In till.' provloiin nootlonn. t Iio t.vpon atui tiia.tn 1 1 tnion of orrorn onoount.oroO 
la porforiiilnf t horiiia 1 - 1 itort 1 a napplti.o. o .1 1 nounnoa . How thono -rrorn tiro 
plaooa In foolofUt oout oxv atU tPoroPy ..vKIo orltorla for ovalnat !n,t t lio on- 
..Iro proooanro. Thnn t iie nolortlon of rio.lollno; aUtorl t Iitn ana a.it a-oo 1 1 oot, I on 
proooauro oan Po oxainino.t In tho I IfP.t of a partloiilar .oooIoo.U- nit nat ion. 

An an 1 11 tint rat I on , annnnio that; wo wish Ho a I no r I rt! nat o Pot worn lltaontono 
ana iiolonIt:o froii; a natolMt.o in a ntiul.v nito tliat. in rolatlvolv flat ;ma whono 
oiitoropn liavo Iron-oxKio nHaInn of approx liaat o ly ono-:a I 1 II aiot m-' t hi oktionn . 


Krror .Ponro< 

iia t.ivl 1 1 1 1* titoannrotaonti 
t.opo.fraphlo oft'oo! n 
pyp'.oal nnrt'aoo-ooat I !if ot'i'ootn 
t.vploal trtinnlont. ofiVotn 


Krrci- Mafnlt; i Mo 

150 Tin 
0 TIU 
.15 TUI 
ICO Tl’J 
:-75 Till . 


Thono ttioantirononv .aiM ttaMo 1 n Imp 1 I ft oat I on Vrrorn aro now oonPInoiI wlt-h 1 ho 
aj:.i!,v:'!n prro.’n. Tho ia>t-al n.vnit. «,*i*t'c'rn t\>t* :*ao!i ,a I f.or I t-Iiia ;irt' : 


A i. . tor It.l ii’i 

i-'onrlor norlon 
.‘•’In ! t o a t f forotioo 
Pa pi a 00 t. rann fortr. 


A p.rrox Inta! o: roi.;i 1 i-;r, 

. ' V50 'I'lU 

•JPO Til' 

550 TUI 


Vvo now know t Intt. wo cannot, a I no r tint nat o Potw.-on f.oolooto mafori a 1 n havltif 
i> _Uto rina U1 no rt. la vil l'••pronoo of lonn I htin V;hl Tli: iinin,: t lio Konrtor-nor Ion t'O ' 
01 ionn t:p.an ‘550 TIU tinlnf t.ho Paplaoo-trar.nfori:; al.torivhn. 

Fl.aiiro li nhown 1 tiorm.a I -I nort: la hi nto.franin ot' varlonn noa 1 imopV - loy ro •!<■•. 

U2. im. Krotii t.ho hlnlofratnn, wo oan aot.ornilno that a t.volcal thermal lnorI.l t 
tor llmo.ntonon In SIMO TI'J aiM for' aol oral t.on in 5‘>00 TIU;' i'.hiin Mio nopar.-it I.on 
in (troat.cr thatrt.ho leant, noonrat.o of tho a 1 .to r 1 1 P.ran . flnco any o!' t.lio th.roo 
analya.ln alporUhmn ctin Pc unoa, tho noleotion of the algorithm will Po av'V.or- 
nliMii Py t.ho cost /i-tYootl vor.osn criteria ana tho Konrlor norlon wn-iia Po 
soloctod. 


In tho hypothotlcal cano. If wo wish ' o a I nor Ira Inato Potwoon naiulntonon 
ana aoloinlt;on, which typically ha\o a thormal Ir.ortla sop.-iratlon of ‘Oil) TUI 
wo mnnt noloot tho i.aplaoo-^t.rannforra a Igor 1 tp.m . " ’ 

.-•i '-'T’ *'■•' V, •'’'■ •■’ Pango of l,tno._nin rooks. 

MU' J on lu't wi'i'U :i i I I noxi roo«< (i‘hvo.l Itv?) nnd n roaiT-o- 

gr.-Unra folnlo rook (gran't.o) In otily .:U0 TIU; tho .•h'l'aration Potwoon rhvniito 
niul a liu'tii mat'lo i*ock (r.abhro) 1? approx Iniatt'ly ^lOO TH). 

tloii a t*liio-i:ral!'.Okl :\>la!c rorU Ti-om a o.oarso-o,ra ! no*-! I'oliUo rov’k can ho 
xKMili'Y...; t!u? Laplaoo-v.ranivt’orm tocluiK^uo, nan.v I V oi'at anti aiiort- 

oiling tho tiiao Intorval. 1’ I nor IpiI nation of a f 1 no-gr.ai nod folnlo rook from a 
ooarno-gralnod in.ifio rook c.'in Po aohlovod with oilhor tho fi iii to-d I f foronoo 
a I ^'.x^r I Ulu c*r tho I.ap)aoo~t rarnTorm a i hm, xiop^Muiinr; rirrt. u|\m} t ho moar- 
uromont. a?ui 'uovio I x-rianax atui aov!o:ui upo:; i-ho Ox^at /orroot i voiior:? . 
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TAB!,F. 2. raramoter valuea used to compare the three algorlUims. - 
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